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Abstract 

Oncolytic viruses are genetically engineered viruses that are designed to kill cancer cells while doing minimal damage to 
normal healthy tissue. After being injected into a tumor, they infect cancer cells, multiply inside them, and when a cancer 
cell is killed they move on to spread and infect other cancer cells. Chondroitinase ABC (Chase-ABC) is a bacterial enzyme 
that can remove a major glioma ECM component, chondroitin sulfate glycosoamino glycans from proteoglycans without 
any deleterious effects in vivo. It has been shown that Chase-ABC treatment is able to promote the spread of the viruses, 
increasing the efficacy of the viral treatment. In this paper we develop a mathematical model to investigate the effect of the 
Chase-ABC on the treatment of glioma by oncolytic viruses (OV). We show that the model's predictions agree with 
experimental results for a spherical glioma. We then use the model to test various treatment options in the heterogeneous 
microenvironment of the brain. The model predicts that separate injections of OV, one into the center of the tumor and 
another outside the tumor will result in better outcome than if the total injection is outside the tumor. In particular, the 
injection of the ECM-degrading enzyme (Chase-ABC) on the periphery of the main tumor core need to be administered in 
an optimal strategy in order to infect and eradicate the infiltrating glioma cells outside the tumor core in addition to 
proliferative cells in the bulk of tumor core. The model also predicts that the size of tumor satellites and distance between 
the primary tumor and multifocal/satellite lesions may be an important factor for the efficacy of the viral therapy with Chase 
treatment. 
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Introduction 

Glioblastoma is one of the most aggressive type of brain cancer 
with the median survival time of approximately one year [1,2] . It is 
characterized by rapid proliferation and high invasion. Glioma are 
resistant to radiotherapy and chemotherapy and eventually recurs 
[3]. Oncolytic viruses (OV) are genetically manipulated viruses 
that can destroy cancer cells but do minimal damage to normal 
healthy tissues [4]. These viruses can preferentially replicate in 
tumor cells, leading to their lytic destruction. Accompanying cell 
lysis the virus particles burst out and proceed to spread and infect 
other cancer cells. While the world's first OV, H101 (Oncorine, an 
OV functionally identical to ONYX-015), was approved by the 
Chinese State Food and Drug Administration, the regulatory 
approval of OVs in the United States and Europe is pending the 
results of randomized and large phase III studies [5] . Despite the 
great potential and regulatory approval of OV for clinical use in 
China, its inefficient dispersal property within the tumor ECM has 
been recognized as a major barrier for its anti-tumor efficacy [5]. 
Tumor ECM plays a pivotal role in inhibiting virus spread [6-1 1] 
and leads to limited viral replication and reduced cytolytic efficacy 
[12]. Major ECM components blocking this viral dispersal in the 



neural ECM include hyaluronic acid (HA) [13,14] and proteogly- 
cans [11,14] and these structural components are also known for 
hindering large therapeutic molecules [11]. 

The use of ECM degradating enzymes could significandy 
increase the efficacy of OV treatment by enabling virus to move 
more freely among uninfected tumor cells. Unlike ECM in normal 
tissue, tumor ECM continuously undergoes remodeling and 
extensive synthesis [15-17]. ECM remodeling is also an important 
part of brain tumor angiogenesis [18,19]. The brain ECM 
contains mainly macromolecules such as glycosaminoglycans 
(GAGs) and proteoglycans (PGs) with a smaller fraction of 
fibrillary glycoproteins such as collagens, fibronectin, or elastin 
[14,19,20]. Tumor ECM inhibits penetration of anti-tumor 
agents, leading to low therapeutic efficacy. High interstitial fluid 
pressure inhibits delivery of the agents by convection through the 
tumor interstitial matrix resulting in a rather passive diffusion as 
the major mode of transport of macromolecules [21,22]. Tumor 
penetration of macromolecules such as IgG is negatively correlated 
with increased collagen content [21,23]. Use of collagenase, an 
enzyme that breaks the peptide bonds in collagen, was shown to 
increase the diffusivity of macromolecules (IgG) [23], to reduce 
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Figure 1 . Overview of ECM degradation in experiments and schematic of the mathematical model. (A) A schematic of ECM degradation 
and Chase ABC treated ECM. (B) In experiments [1 1], viral spread was evaluated 4 days after OV treatment through detection of the GFP-positive 
infected cells. Images in upper and lower panels show 3 representative spheroids (n = 6/group) treated with vehicle (Chase-) and Chase-ABC (Chase+), 
respectively (Figure 1D in [11]). Chase-treated viruses spread more efficiently throughout the spheroid compared to control-treated spheroids. 
Glioma U87AEGFR spheroids were cultured on organotypic brain slices for 72 hours followed by treatment with purified Chase-ABC or vehicle for 
24 hours before or after infection with rQNestin34.5. (with permission from the following article: Dmitrieva N, Yu L, Viapiano M, Cripe TP, Chiocca EA, 
Glorioso JC, Kaur B., Chondroitinase ABC l-mediated enhancement of oncolytic virus spread and antitumor efficacy, Clin Cancer Res, 17(6):l362-72 (2011). 
(C) A schematic of the mathematical model. Three types of glioma cells (uninfected (a), infected (v), and dead («) cells) and ECM component (p) 
consist of a growing tumor mass with a free boundary T(t). Viruses (v) from outside (or inside) the tumor mass penetrates through the network of 
ECM components within the the tumor mass to infect glioma cells (a). Diffused Chase-ABC (C) degrades ECM components and help virus infiltrate the 
narrow intercellular space between cells within the tumor mass. 
doi:1 0.1 371/journal.pone.01 02499.g001 



interstitial fluid pressure, and to enhance the transvascular 
transport through convection [24]. 

Protease- or hyaluronidase-mediated digestion of the ECM can 
enhance intratumoral penetration [25], and MMP-enhancing 
oncolytic relaxin is known to increase anti-tumor efficacy of virus 
spread [7] . Relaxin is a peptide hormone that is able to reduce the 
synthesis and secretion of interstitial collagens and increase the 
level of metalloproteinase (MMP) and oncolytic adenovirus 
expressing relaxin was shown to promote OV dispersal through 
the tumor ECM [7] . However, ECM degrading enzymes need to 
be used with caution since intracranical use of degrading enzymes 
can cause serious complications. For example, hemorrhagic 
necrosis of brain can be induced by collagenase-mediated ECM 
disruption, brain proteases are associated with neurodegenerative 
diseases [26], and hyaluronidase can promote optic glioma growth 
through astrocytic reactivity. [27]. 

Secreted and membrane-bound chondroitin sulfate proteogly- 
cans (CSPG) linked to extracellular hyaluronan is one of major 
components of the ECM in the brain [13]. CSPGs include 
versican, brevican, phosphacan, and NG2. CSPGs also play a role 
as axon growth inhibitory molecules that are present in the glial 
scar, and are responsible in part for regeneration failure of scar 
after damage to the CNS [28]. In CNS tumors, increased 
expression of CSPGs is associated with tumor growth, invasion, 
and angiogenesis [29], and is responsible for diffusion-limiting 
properties [1 1]. Accumulated CSPGs in the glioma ECM increase 



tortuosity of the extracellular space and interstitial pressure within 
the tumor [11], leading to poor transport of large therapeutic 
agents [22]. Hence enzymatic manipulation, lowering the CSPG 
levels, would be a way of improving interstitial transport of 
therapeutic agents into the tumor. In this connection, ECM 
macromolecules have been considered as potential therapeutic 
targets for adjuvant therapy [30,31]. For example, tumor 
proliferation and dispersion were inhibited when these molecules 
were reduced or inhibited by blocking antibodies against versican 
[30] or interfering RNA against phosphacan [32]. CSPGs can also 
be degraded by MMP-1 and MMP-8, leading an increase in 
hydraulic conductivity and particle diffusion in solid tumors [33]. 

Chondroitinase ABC I (Chase-ABC) is a bacterial enzyme that 
can remove Chondroitin sulfate glycosoamino glycans from 
proteoglycans without any deleterious effects in vivo [1 1] and 
has been studied for its effect on neuronal regeneration after 
injury. Chase-ABC has been widely used to enhance regeneration 
of injured axonal tracts due to its a long-lasting 'loosening' effect 
on the ECM scaffold [28,34]. Based on extensive positive evidence 
from preclinical models, Chase-ABC I have recendy been used in 
Phase I/II trials for treatment of patients in Japan [35] . In an 
effort to investigate the use of bacterial chondroitinase to enhance 
anticancer therapy, Dmitrieva et al. [1 1] examined the effect of 
Chase-ABC on tumor ECM, OV spread, and efficacy. In the 
experiments, three-dimensional glioma spheroids placed on 
cultured brain slices were utilized to evaluate OV spread; see 
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Figure 1A. Replication-conditional OV expressing Chase ABC 
(OV-Chase) was tested for spread and efficacy in glioma spheroids 
and glioma xenografts were utilized to compare anti-tumor 
efficacy of OV-Chase, HSVO_ (control) and PBS. Dmitrieva et 
al. [1 1] found that OV spread in glioma spheroids grown on brain 
slices were significandy enhanced by Chase-ABC treatment, that 
cell migration or invasion were not enhanced by OV Chase 
treatment. Together these results suggest that degradation of ECM 
by OV expressed bacterial Chase-ABC was the key factor in 
enhancing OV spread and anti tumor efficacy. It has been shown 
that virus spreads more efficiently in spheroids compared to 
control-treated spheroids when they are treated with Chase [1 1]; 
see Figure IB. 

Kim et al. [36] recently developed a mathematical model of 
invasive glioma cells in 3D tumor spheroids. The model was 
shown to be capable of reproducing migration patterns of glioma 
cells in in vitro experiments, exhibiting, in particular, dispersion 
and branching of cells. The model included MMP activity and 
glucose levels as well as chemotaxis, haptotaxis and cell-cell 
adhesion forces. The rapid migration of cells is caused primarily by 
the chemotaxis forces including glucose levels. There are several 
mathematical models of brain tumors based on reaction-diffusion 
process [37-41], cellular automata [42], the transition between 
migration and proliferation using kinetic or diffusion models [43— 
45], and intracellular signaling pathways [1,46,47]. A general 
review on mathematical models of tumor growth can be found in 
[48,49]. 

In this paper we develop a mathematical model based on the 
schematic diagram of Figure 1C. The model consists of a coupled 
system of partial differential equations involving the following 
variables: 



x= density of uninfected cancer cells, 
y = density of inf ected cancer cells, 
n = density of dead cancer cells, 
E= concentration of tumoral chondroitin sulfate proteoglycans (CSPG), 
v= density of virus particles, 
C= concentration of chondroitinase ABC (Chase ABC), 
p = concentration of tumor ECM material. 



It will be shown that the model simulations for a spherical 
glioma are in qualitative agreement with experimental results of 
[1 1]. It is a common practice to inject the virus into serval different 
locations within the tumor. The mathematical model can be used 
to determine preferred locations for such injections. In Materials 
and Methods section, we present the mathematical model, and in 
Results Section we simulate the model, and then also perform 
several experiments and propose some hypotheses. Parameter 
estimation, nondimensionalization of the model and numerical 
scheme of dimensionless governing equations are given in Text S 1 . 



growth; in Section 2.6, we will derive an equation for u. The 
equations for x,y and n are taken from the model by Friedman et 
al. [50]; 
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where u is the velocity, X is the proliferation rate of uninfected 
tumor cells, /? is the infection rate, & is the infected cell lysis rate, 
and fi is the removal rate of dead cells. We note that treatment by 
Chase-ABC does not affect these equations. We also note that 
Friedman et al. [50] included the effect of the immune system and 
immune suppressive treatment on the efficacy of the OV injection; 
in the present paper we ignore this aspect of the treatment, 
concentrating only on the efficacy derived from Chase-ABC 
treatment. 

Concentration of tumoral chondroitin sulfate 
proteoglycans (CSPG) (E) 

Secreted and membrane-bound chondroitin sulfate proteogly- 
cans (CSPG), described below in Section 0, linked to extracellular 
hyaluronan form a major component of the ECM in the brain 
[13]. It was found that CSPGs detected in glioma cell lines and 
tumor samples interfere with OV spread in glioma spheroids [1 1]. 
The ability of Chondroitinase ABC (C) to degrade the tumoral 
chondroitin sulfate proteoglycans (CSPG) is the key regulation 
factor in this paper. As in [11], we hypothesize that Chase-ABC- 
mediated digestion of glioma CS-GAGs degrades glioma ECM 
and opens up space for viruses to diffuse faster. 

The conservation of mass for the CSPG concentration gives 



3E 

~dt 



= -V-J e + Pe, 



(4) 



where Je is the flux for mobility and Pe is the CSPG production/ 
death rate. The CSPG is convected by the fluid velocity u, 



Je = Eu. 



(5) 



Materials and Methods 

In this section we develop a system of partial differential 
equations with initial and boundary conditions, for the variables 
introduced in Section 1, in a spherical tumor Q(t) and boundary 

r(f). 

Cell density (x,y,ri) 

In order to describe the time evolution of densities of different 
types of tumor cells (uninfected, infected, and dead cells), we take 
into account cell proliferation of tumor cells and the virus infection 
of tumor cells, followed by the clearance process, and the passive 
movement due to the velocity field u generated from tumor 



CSPG degradation rate is given by 

P E =f(E,C) (6) 

where f(E,Q describes the degradation of the ECM (CSPG) by 
Chondroitinase ABC (C(x,t)). 

Chase-ABC degrades the key ECM component CSPG(£). We 
assume that the ECM degradation rate is proportional to both the 
ECM concentration and Chase concentration, 
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f(E,Q = 



1 17 cm 

/L4] Hi C- 

K'" + C" 



(7) 



where M\ is a degradation rate of the ECM component, CSPG. 

In summary, we have the following governing equation for the 
CSPG concentration 



d_E 

Tt 



+ V(Eu)-- 



A41E C m 



Removal by Chase-ABC 



(8) 



This equation describes the time evolution of the CSPG density 
after Chase-ABC-assisted degradation of CSPG with the passive 
movement from the velocity u. 



where y is the clearance rate of viruses, and I a (') is the indicator 
(characteristic) function over a domain A. 

In summary, we have the following governing equation for the 
virus density, 



dv 
8i' 



v-( 



K E +E% 



Random motility 
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(14) 



Chondroitinase ABC (Chase-ABC; C) 

We next derive the time evolution of Chase-ABC concentration 
in response to injection, secretion and ECM degradation, and 
natural decay, in the system. From conservation of mass for the 
concentration of Chase-ABC (=C), we get 



Density of free virus particles (v) 

In this section we describe the movement and replication/decay 
process of virus particles within the complex ECM structure, 
CSPG, in the microenvironment. Treatment with Chase-ABC 
improves viral spread [11]. We assume the virus can spread in 
CSPG concentration-dependent manner and duplicate and 
undergoes apoptosis. From conservation of mass for density of 
free virus particles ( = v(r,t)), we get 



dv 
It 



= -V-J v + P v , 



(9) 



where J v is the flux for mobility and P v is the birth/death rate. 
The only contribution to the flux J v is flux from random motility 
of the free viruses: 



J r =-D(E)Vv 



(10) 



where D(E) is the random motility that depends on ECM (CSPG) 
concentration E and quantifies how virus particles outside infected 
cells can freely move around in less dense ECM (CSPG). We use 
the following phenomenological form 



D(E)- 



Ke+E° 



-D 



PBS 



(11) 



where Dpss is the diffusion coefficient in the PBS (reference) case, 
Ke is an inhibition parameter, A51 is a scaling factor, and E s aw is 
the average local density of the ECM concentration over a local 
neighborhood, i.e., 



8C 



= -S/-J c + Pc 



(15) 



where Jc is the flux and Pq is the production/ decay rate of the 
Chase-ABC. The only contribution to the flux Jc is flux from the 
diffusion process: 



J C =-D C VC 



(16) 



where Dc is the diffusion coefficient of the Chase-ABC. The 
Chase-ABC is applied in the far away field and degrades the key 
ECM component CSPG (E). Chase-ABC may be secreted by 
infected tumor cells through virus [5 1] or be directly applied to the 
field surrounding a growing tumor in addition to the injection of 
virus [1 1]. Taking into account the natural decay of the molecule, 
we set 



P c =g(v,y,Sl)-X a EC-l 62 C 



(17) 



where g(v,y,Q) is a function that describes a production of 
Chondroitinase, which depends on virus, infected tumor cells, and 
possibly other factors in the tumor Q(t), ^3 is the ECM 
degradation parameter, and A42 is the decay rate. 

In summary, the governing equation for the Chase-ABC is as 
follows: 



= V (Z)cVQ + g(v,y,Sl) - laEC 

Diffusion Secretion Degradation 



^6iC . (18) 

Natural decay 



1 



|£,(x) 



E(z,t)dz, 



(12) 



where B s (x) = {zeR 3 : \x-z\<s} is the ball in 3D (or a disc in 
2D) with center x and a. sensing radius s, and |-B s (x)| is the volume 
of -S.s(x). Here we use the convention that E(z,t) = 0 if z is outside 

A 

the tumor Sl(t). When E = 0, this gives D(E)= —D PBS . 

Ke 

We assume free virus particle can replicate with replication rate 
b in the presence of infected cells (y) with lysis (or death) rate <5 
within the growing tumor [50] 



P y = bSyI m -yv 



Tumor ECM material (p) 

Glioma cells destabilize the normal cell-matrix structure by 
overexpressing or generating their own altered neural ECM 
molecules [52,53], creating a disorganized ECM structure that 
promotes tumor-cell movement and proliferation [13]. We assume 
that only proliferating uninfected tumor cells (x) secrete tumor 
ECM materials that undergo remodeling. By conservation of mass, 
we get 



dp x m 7- p 



(19) 



(13) where X42 is secretion rate of ECM material by uninfected tumor 
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cells, OT2 is the Hill-type function exponent, K x is the Hill-type 
function parameter, p 0 is the maximum capacity of the tumor 
ECM material. 

Growing spheroid 

The tumor spheroid either grows or shrinks in response to 
growing mass from proliferative cells and added ECM materials, 
and decreased mass from the ECM degradation process. We 
model this process by the Stokes equation with added material 
acting as added pressure: 

~Au + Vp= (2Q) 

where p is the pressure and r\ is the viscosity of fluid. The first two 
terms, third, and fourth terms on the right-hand side of equation 
(20) represent additional total change in the density of cells, new 
material of the ECM, and degraded ECM (CSPG) components, 
respectively. By conservation of mass, div u is equal to the 
expression in parenthesis on the right-hand side of equation (20). 

Boundary Condition. The governing equations (l)-(20) is 
coupled with the following boundary conditions 

^ = «, on r(0 (21) 



^=0, onaC2 0 , (22) 



— =0, on <5C2 0 . (23) 
dv 

where v is the outer normal vector and X = X(t) represents the 
free boundary. Parameter and reference values are given in 
Tables 1 and 2. 

Results 

In next two sections we take g(v,y,Q.) = Q. 

Results and comparison with experimental data 

We first compare our simulation results with experimental data 
in [11]. In experiments [11], various human glioma cell lines 
(U343, U87, U87DEGFR, LN229, Gli36DEGFR-H2B-RFP, 
XI 2) and Vero cells were used in in vitro and in vivo expriments. 
OV-Chase is an oncolytic HSV-1 -expressing bacterial Chase-ABC 
driven by a promoter within the backbone of rHsvQ_ (control). In 
vitro studies, in order to see the spread of rHsvQ, and OV-Chase 
in 3-dimensional glioma spheroids on brain slices, the spheroids 
were infected with 10 4 pfu (plaque forming units) of rHsvQ, or 
OV-Chase. Then, infected areas for the cases of Chase(— ) and 
Chase(+) were calculated by both visualization and quantification 
of infected cells by appearance of OV-encoded GFP over time. In 
in vivo studies, the glioma cells implanted stereotactically at a 
position 2 mm lateral to bregma at a depth of 3 mm for animal 
studies and animals were inoculated with control (PBS) or virus in 
the absence or presence of Chase-ABC ten days after tumor cell 
implantation. Intratumoral injections of HBSS or the indicated 
OV was performed for the mice bearing subcutaneous tumors 



with the initial average volume of 80 to 150 mm 3 . Seventy-two 
hours after viral infection, the amount of infectious viral particles 
and infected areas were quantified by a standard method. 
Figures 2A-2L show spatial profiles of uninfected cells (x), infected 
cells 0), dead cells («), virus (v) and CSPG ECM (E) at t = 24, 60h 
without and with Chase treatment. Without Chase treatment, 
CSPG ECM degradation is inhibited (Figures 2Q,2R) and virus 
spread (Figures 2M,2N) is limited due to thick ECM components. 
Figures 2M,2N show localized virus in the annulus subdomain. 
Infected tumor cells and dead cells are localized at the periphery of 
the spheroid due to this inhibited virus spread toward the center of 
the spheroid. There are large portion of uninfected tumor cells at 
the center of the tumor spheroid. On the other hand, Chase 
treatment promoted CSPG ECM degradation and decomposition 
(Figures 2S,2T) and led to fast and effective spread of virus 
(Figures 20,2P) and more infection of tumor cells near the 
spheroid center (Figures 2G,2H). One also observes significant 
reduction in uninfected tumor cell density (Figures 2C,2D) and 
widespread distribution of dead tumor cells at the center of the 
tumor spheroid (Figures 2K,2L). 

In order to compare the computational results with experimen- 
tal data, we calculated the relative infected area using density 
distribution of infected tumor cells i.e., we use the infected area 
relative to the total area of the glioma spheroid. Figure 2U shows 
relative infected areas with (black) and without OV-Chase 
treatment. This result is in good qualitative agreement with 
experimental data in [1 1] shown in Figure 2V. The experiments 
in [11] and the simulations in this paper demonstrate that OV- 
chase treatment improve the infection of virus within the tumor 
thus decreasing the burden of cancer mass. 

The role of the microenvironment 

The microenvironment plays a significant role in regulating the 
spread and treatment efficacy of many types of cancer [36,54,55] 
including glioblastoma [3]. Chase diffusivity which clearly affects 
the tumor microenvironment may depend on many factors such as 
composition of ECM materials in in vivo system. In Figure 3, we 
investigated effect of diffusivity of Chase on virus spread. The 
relative infected area is increased when diffusivity of Chase is 
increased. Almost all tumor cells at the center of tumor mass were 
infected when diffusivity is 1.5-fold larger (1.5Z>c) than was in the 
control case (Dc)- On the other hand, when diffusivity of Chase 
was decreased, the virus spread was severely impaired. For 
example, the infected area was decreased by 50% and 37% and 
only the cells on the periphery of the tumor mass were infected 
when the diffusion coefficient of Chase (Dc) is 10- and 100-fold 
smaller, respectively. This implies that virus spread would be 
significantly less effective when a tumor is located in a harsh 
microenvironment where ECM component may be densely 
packed or quickly remodeled. 

Figure 4 shows virus spread at t = 24,48, 60A when a combina- 
tion of different diffusion coefficients on left and right half of the 
square domain (second column in Figure 4). Patterns of Infected 
areas for different combinations show that virus spread is 
significantly inhibited in the area of low diffusivity. When 
diffusivity of Chase is fixed on the left half of the domain but 
decrease in the right half of the domain, most tumor cells on the 
left side were infected but not many cells near the center of the 
tumor mass on the right hand side were infected. For example, in 
Figure 4C, base value of diffusivity was assigned on the left half 
(Dc = Dc) an d 100-fold smaller values of diffusivity of Chase was 
assigned on the right half [Dc = D c/100). This difference in 
diffusivity leads to different patterns of virus infection on the left 
and right half of the domain. While virus infection is completed 
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Table 1. Model parameters. 




Parameter 


Description 


Dimensional value 




Ref 


Diffusion coefficients/Random motility 




Random motility of virus (PBS) 


3.89 x 10 


- 2 {mm 2 jh) 




[89], tw 


D c 


Diffusion coefficient of Chase 


3.89 x 10 


- 4 {mm 1 jh) 




tw 




virus random motility controlling parameter 


2.25 x 10 


- mm /g 




tw 


K E 


Inhibition parameter of virus motility from ECM 


l.Ox 10" 


°mm /g 




tw 


Production/remodeling rates 


A 


Proliferation rate of tumor cells 


1.54 x 10 


(1/A) 




[50], tw 


Xt) 


carrying capacity of uninfected tumor cells 


-X* 






tw 


(< 


Infection rate 


7.0 x 10" 


10 [mm jh virus) 




[50] 


b 


Burst size of infected cells 


50 (virus/ cell) 




[50] 


Ml 


ECM production rate from uninfected tumor cells 


6.0 x 10- 






[36], tw 


K x 


Hill-type coefficient in tumor ECM 
remodeling/reconstruction by uninfected tumor cells 


=x* 






tw 


Po 


carrying capacity of tumor ECM 


= P * 






[36] 


1112 


Hill-type parameter in ECM remodeling/reconstruction 
by uninfected tumor cells 


it 






tw 


AC 


Secretion rate of ECM degrading enzyme 


i.OxW'mU/Qi.g) 




tw 


Inhibition/degradation/decay rates 


y 


Clearance rate of viruses 


1.8 x 10- 


3 tt/h) 




[50] 


s 


Infected cell lysis rate 


1.2 x 10" 


2 a/h) 




[50], tw 


fi 


Removal rate of dead cells 


1.2x10- 


2 a/h) 




[50], tw 




ECM degradation rate 


9.0 x 10 1 


a/m 




tw 


mi 


Hill-type parameter in ECM degradation by Chase 


it 






tw 


K c 


Hill-type coefficient in ECM degradation by Chase 


= C* 






tw 




Chase reaction rate due ECM degradation 


6.0 x W i mm 3 /(h.g) 




tw 




Natural decay rate of Chase 


6.0 x 10- 


3 fX/h) 




[28], tw 


Units are indicated in parenthesis ( ) in the third column. 
1 =dimensionless value. 
*tw = estimated in this work. 
doi:10.1371/journal.pone.0102499.t001 










Table 2. Reference value used in the model. 












Var Description 


Dimensional value 






Ref 




L Length scale 


3 mm 






[11] 




T Time scale 


1 .67 h 






tw 




D Characteristic diffusion coefficient 1.5 x \0~ 5 cm z /s 






tw 




x* Uninfected cell density 


\0 6 cells/mm 3 






[50,92] 




y* Infected cell density 


= X* 






[50,92] 




n* Dead cell density 


= n* 






[50,92] 




v* virus concentration 


2.2 x 10 8 v/™.y/ mm 3 






[50,92] 




p* Concentration of other ECM materials l.Omg/cm 3 






[100,101] 


E* ECM (CSPG) concentration \.0mg/cm 2 






[94-100] 


C* Concentration of Chase-ABC C* = SOmU/ml 






[28,34,93] 


tw = estimated in this work. 
doi:10.1371/journal.pone.0102499.t002 
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W/O Chase ABC W/ Chase ABC Variables 

t=24h t=60h t=24h t=60h 




(M) (n) (o) (p) n 2 




Figure 2. Dynamics of the model and comparison with experimental data. (A-T) Dynamics of the model in the control case (without chase 
treatment) and with chase treatment. Profiles of densities of uninfected, infected, and dead cells, virus, and CSPG ECM, are shown at 24, and 60 hours 
postinfection. Without chase treatment, virus spread is limited due to thick ECM components. With Chase treatment, ECM degradation is improved, 
leading to enhanced virus spread. The enclosing box size is [0, 1] x [0, 1]. ac = 0. (U, V) OV-Chase has improved virus spread in glioma spheroids. 
Relative infected areas [i.e., up to scale) at t = 24,48, 60/i were calculated from model simulations (U) and experiments (V). In simulations, viral spread 
was calculated as the infected area relative to the total area of the glioma spheroid. The black and white bars represent the results with and without 
ChaseABC, respectively. ).c — 0. In experiments, viral spread was calculated as the GFP-positive area relative to the total area of the spheroids in [11]. 
doi:1 0.1 371 /journal.pone.01 02499.g002 



throughout the left half, only tumor cells on the periphery of the 
spheroid was infected in the right half of the domain. Similar 
patterns were observed in the case of 10-fold smaller Dc value on 
the right in Figure 4D (D C = D C on the left; D C =D C /10 on the 
right). When the diffusivity is decreased on the left as well and 100- 
fold smaller diffusivity was fixed on the right (Dc = Dc/ 100), the 
infected area on the left also is decreased (Figures 4A-4B) as well. 
Figures 4E-4H show relative infected areas on the left, right, and 
whole domain, respectively, for the four cases in Figures 4A-4D. 
The simulation results in Figure 4 have implications for the 
efficacy of virus spread in the brain when chase components are 
injected a site near a tumor since brain tissue including white and 
gray matter is heterogeneous in nature. 



Testing hypothesis for real tumors in vivo and optimal 
solution of OV treatment 

Muir et al. [51] were able to induce secretion of bacterial 
chondroitinase ABC from mammalian cells by modifying N- 
glycosylation sites. In this section, we assumed that infected tumor 
cells secrete the enzyme since the virus use the cell body to create 
this enzyme inside the infected cell. In this section, we use our 
model to test the results of treatment of tumor in vivo by 
genetically engineered virus that secrete ECM degradating 
enzyme through infected cells instead of injecting enzyme C, 
which acts passively outside the growing spheroid. In order to take 
into account the effect of secretion of enzyme by virus through 
infected cells we set in the equation (18) 
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Figure 3. Effect of Chase diffusion on virus infection. (A-E) Viral spread were shown at 60 hours postinfection for various values of Dc 
(D c = Z>c/100 in (A), Dc/10 in (B), D c /2 in (C), D c in (D), 1.5D C in (E)). The default value of D c , D c , is 7.2 x 10~ 5 . The enclosing box size is 
[0.15,0.85] x [0.15,0.85]. (F) Relative infected area for each case in (A-E). 
doi:1 0.1 371 /journal.pone.01 02499.g003 



g(v,y,Q) = X c y 



(24) 



where Xq is the secretion rate of the enzyme from infected cells. 

Figure 5 shows profiles of infected tumor cells, at time t = 60/?, 
in response to various secretion rate of the enzyme by infected cells 
(Ac = 0, 0.005, 0.01, 0.02, 0.03). An increase in indirect secretion 
of the enzyme by infected cells effectively promotes penetration of 
the enzyme toward the tumor center and degradation of the ECM 
components: as Xc is increased, the infected area is increased. For 
example, while tumor cells at the center of tumor mass were not 
infected at t = 60/j in the the control case, almost all tumor cells 
were infected even at the tumor center when Ac = 0.03. 

We next consider a situation where be virus injection takes place 
at different locations. Figures 6A-6C show the efficacy of 
redistributed injection for the OV treatment with/without 
Chase-ABC relative the control case. We first inject the 15% of 
total OVs in the center of the growing glioma with/without 
treatment of Chase and inoculate the remaining OVs (85%) on the 
periphery of the tumor along with Chase treatment. For control 
case, same amount of total virus was injected outside the tumor. 
Profiles of infected cells at t — 9, 18, 23, 60/j for control case (A), 
case without chase (B) and case with chase (C) were shown in 
Figures 6A-6C. While virus penetration begins on the periphery of 
the spheroid and stops in middle-way toward the center at final 
time (t = 60h) in the control case (Figure 6A), viral infection at the 
center of the spheroid spread outward aggressively in the case of 
viral injection at the center (Figures 6B-6C). In the case of Chase 
treatment at the injection site at the center of the spheroid 
(Figure 6C), viral infection is even more aggressive and almost all 
tumor cells were infected at final time (t = 60h). Figure 6D shows 
time courses of total relative areas for the case of control (black 



solid; Figure 6A), without Chase (blue dotted; Figure 6B), and 
with Chase (red dotted; Figure 6C). Injection strategies on 
different locations (red,blue) improves the tumor infection 
compared to the control (black solid). Treatment of Chase at the 
tumor center (red dotted) leads to even faster viral spread. 
Figure 6E shows time courses of relative areas of outer ring (solid) 
and inner region (dotted) for the case without Chase (blue, 
Figure 6B), and with Chase (red, Figure 6C). While infection areas 
of outer ring for both cases with/without Chase treatment are 
same at / = 24, 48, 60, 8O/1 (blue solid, red solid), the case with 
Chase treatment (red dotted) at the center of tumor spheroid leads 
to faster spread of virus compared to the case without Chase (blue 
dotted). The injection amount of virus is limited due to safety 
concern to the patients from increased pressure inside the brain. 
Therefore, given MRI images of patient brain, we would like to 
optimize the efficacy of the treatment with limited number of virus 
per total injection, by carefully choosing injection sites. Our 
simulation results (Figure 6) which show improved efficacy by viral 
injections at both the center of the tumor and on its periphery, are 
somewhat supported by recent clinical data reported in Markert et 
al. [56]. Indeed, Markert et al. [56] investigated the efficacy of OV 
treatment in six human patients by two injections of OV. The first 
inoculation of 13% of 1.15 x 10 9 pfu of G207 in the center of the 
tumor mass was followed by a surgical resection of the tumor mass. 
The remaining virus (87%) was injected into brain tissue near the 
boundary of resected tumor area. In this study, virus was injected 
in tumor tissue by a catheter that was left in place, and the tumor 
resected en bloc 2-5 days later. Tissue was carefully divided into 
three sections a distal slice, a more proximal and the most 
proximal referred. Evaluation of HSV polymerase DNA in the 
tissue showed a decrease of 1 to two logs between the distal and 
proximal tissue in four of the six patients. 
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Figure 4. The role of the microenvironment in regulating the virus spread. (A-D) Virus spreads were shown at 24, 48, and 60 hours 
postinfection when a tumor was imposed in the different microenvironment. Slower transport of virus in a harsh microenvironment (slower diffusion) 
decreases virus spread efficacy. Different diffusion coefficients of Dc were prescribed on the left ([0,1] x [0,0.5]) and right ([0.5,1] x [0,1]) half of the 
„D C D c . , . , , , . . fbc b c , ^ . .be be 
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on the left and right in (C); Py^.^jTr) on the left 



domain [0,1] 2 Kw.jw) on the left and ri 9 ht in ( A >; ^'ifjrj' on the left and right in ' B ' ; \ 
and right in (D);). The enclosing box size is [0,1] x [0,1]. The relative infected areas at 60 hours of (A), (B), (C), and (D) are 0.4397, 0.5845, 0.6743, and 
0.7258, respectively. (E-H) Dynamics of virus spread in different microenvironment in (A-D). Relative infected areas on the left half (blue, dotted), right 
(red, dotted) half, and whole (black, solid) tumor spheroid for four cases in Figures 4A-4D ([Left, right] = [b c /W,b c / W0] in (A), [Z> C /2,Z> C /100] in (B), 
[bc,bc/W0] in (C), [b c ,b c /W] in (D)) 
doi:1 0.1 371 /journal.pone.01 02499.g004 



Testing hypothesis on efficacy of a viral therapy on 
eradicating the invasive glioma cells 

High invasiveness of glioma cells is the major challenge in 
treatment of glioma due to regrowth of the tumor after 
conventional treatment options such as surgery and chemo- and 
radio-therapy. Innovative therapeutic approaches of targeting 
these invasive cells are necessary in order to improve clinical 
outcome [57]. Here we test some hypotheses on the efficacy of a 
possible viral therapy on eradicating invasive glioma cells that shed 
from the main tumor core. For this purpose, we assume that the 
Chase-ABC is delivered by injection at different sites (Q»), so that 



Total Flux of Chase-ABC = V h in I aB _ 

:_1 J 



(25) 



7=1 



where iV,„ is the number of injection sites and hj n is the injection 
strength at £2#.. For simplicity, we assume an initial configuration 

of a spherical tumor in the brain tissue with invasive cells at 
random locations outside of the tumor core as in Figure 7A. The 
injection sites were assumed to be located on the periphery of the 
spherical tumor. Figures 7B-7H show the density of the infected 
tumor cells for injections of Chase-ABC at different sites, with the 
same strength at each site after r = 60A. In Figure 7B the Chase- 
ABC is injected near the north east (NE) corner, resulting in more 
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Figure 5. Effect of X c on virus infection. (A-E) The profiles of infected cells for various X c (X c = 0 in (A), 0.0005 in (B), 0.01 in (C), 0.02 in (D), 0.03 in 
(E)) were shown at 60 hours postinfection. The enclosing box size is [0.15,0.85] x [0.15,0.85]. (F) Infected area at t = 60h for various Xc in (A-E). As Xc is 
increased, the infected area is increased. 
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virus-infected area (i.e., more virus-infected tumor cells) in the NE 
area than elsewhere. On the other hand, since the extracellular 
virus do not move freely in the south west (SW) corner (where no 
Chase-ABC was injected), the virus are more localized, resulting in 
the higher concentration of infected tumor cells, but in very 
narrow strips. Figures 7B-7H show that by strategically injecting 
increasing number of Chase-ABC, the number of infected cells 
increases, not only inside the tumor but also among the invasive 
cells We note however that very little improvement occurs after 
injections at more than four sites. This is so because the Chase- 
ABC at the four sites manages to diffuse quite effectively in the 
entire region after 60 hours, and, furthermore, the effectiveness of 
Chase-ABC in the removal of CSPG reaches saturation level Am 
as the concentration of Chase-ABC increases to infinity (see 
equation (8)). We conclude that for a tumor as in Figure 7A, it is 
not necessary to inject Chase-ABC in more than four strategically 
chosen sites. Larger tumors, or tumors with more invasive cells 
distribution may require additional injection sites. The relative 
infected areas for invasive (white), core (gray), and total glioma 
cells in seven cases are shown in Figure 8B. As the number of 
injection sites is increased, the infected population of both invasive 
and proliferating tumor cells inside the core is increased overall. 
The model predicts high improvement on the efficacy of the 
treatment when Ni„ is increased from 1 to 4 (white bar in 
Figure 8B). However, the efficacy of the treatment on infecting the 
infiltrating glioma cells is not significantly improved if Nj„ is 
further increased (Ni„ > 4) due to the saturation phenomenon 
mentioned above. Figure 8A shows the relative infected area of 
both infiltrating and proliferative tumor cells for various number of 
injection sites (Njn) and injection strength (/z,-„). The efficacy is 
increased as N/ n is increased for fixed /z,„ and as h m is increased for 
fixed N in . However, because of saturation effect, there is no 
significant increase in the total infected area if A,„ > 4 or if hi n is 
increased beyond a certain limit. Since too many injections or high 



concentrations of Chase-ABC may harm normal cells, our model 
suggests that the choice of a small number of strategically located 
sites (e.g. four) with limited concentration of Chase-ABC is the 
optimal strategy that would be effective in eradicating tumor cells 
not only in the visible core tumor but also in the isolated cluster as 
long as the cluster of invasive cells are not too far from the main 
tumor bulk. 

Virus could also be transported to brain tumors through the 
vascular route [58]. We test the efficacy of Chase-ABC on killing 
both infiltrating cells and tumor bulk when virus supply is through 
intravenous injections at different sites (Qb>). Thus we have the 
following virus supply via the intravenous injection: 



Total Flux of virus = S~] h\ n I a 



(26) 



where is the number of blood sites for virus transport and h v jn is 

the strength (or concentration) of the viruses at i2gv. Invasive 

,/ 

tumor cells were initially distributed at random locations outside 
the tumor bulk (lower panel in Figures 9A,9D). In the absence of 
Chase-ABC, viral infection are not effective due to limited virus 
spread for both invasive cells and proliferating cells in the main 
tumor bulk when viruses were supplied through one (Figure 9B) 
and five (Figure 9C) blood vessel sites. The infection rate for 
infiltrating cells is increased as NJ n is increased (Figure 9B — > 
Figure 9C). However, there is no significant improvement in 
overall penetration of the main tumor bulk due to lack of Chase- 
ABC. On the other hand, there is a significant improvement in 
infecting the bulk tumor when Chase-ABC was injected at four 
sites (Figures 9E,9F). Upper panels in Figures 9B,9C,9E,9F show 
the initial distribution of viruses. In Figure 9B a higher viral 
infection is also observed in the upper part of the domain in the 
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Figure 6. Effect of intratumoral injections of virus. Partial portion of virus (%1 5 of virus) was injected at the tumor center with (B)/without (C) 
Chase treatment and remaining %85 of virus was injected outside the tumor with Chase treatment at the tumor center For control case, same 
amount of total virus was injected outside the tumor. (A-C) Profiles of infected cells at t = 24,48, 60, 80/i for control case (A), case without chase (B) and 
case with chase (C). Dc = 1-2 x 10 ~ 6 while other parameters are fixed as in Table 1. (D) Time course of total relative areas for the case of control (black 
solid; Figure 6A), without Chase (blue dotted; Figure 6B), and with Chase (red dotted; Figure 6C). Injection strategies on different locations improves 
the tumor infection relative the control. With treatment of Chase at the tumor center, relative infected area is even further improved. (E) Time course 
of relative areas of outer ring (solid) and inner region (dotted) for the case without Chase (blue, Figure 6B), and with Chase (red, Figure 6C). 
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presence of Chase-ABC due to the location of the virus injection 
site. Thus, presence of more blood vessels close to the invasive 
glioma cells and bulk tumor naturally increases the efficacy of 
killing infiltrating tumor cells by viral infection (Figure 9E— > 
Figure 9F). As in the previous test above, the larger number of 
blood sites as a source of viral injection (NJ H ) does not significantly 
improve the efficacy of viral spread and infection for both invasive 
cells and the bulk tumor (not shown here). Too large numbers of 
viruses may harm normal cells and there exist a limit for virus 
injection. Especially in the case of intravenous injections, control 
over how much viruses can be transported is rather limited due to 
the different microenvironment that the tumor resides in. In brain, 
entry of even small biological agents into brain tissue is usually 
more complex than other organs because of the blood-brain 
barrier (BBB) [59,60]. Despite successful transport of virus via 
blood routes and viral infection of gliomas (VSVrp30a; [58]), the 
delivery of viruses may depend on many factors such as 



biochemical conditions of blood vessels, tissue composition and 
geometry. Therefore, given blood distributions of the patients, one 
has to use optimized patient-specific injection strategies with 
minimal side effects, i.e., without killing normal healthy brain cells, 
in order to eradicate both invasive cells and the tumor bulk. 

Effect of virus injection with Chase-ABC on killing primary 
and secondary tumor cells 

The experimental results reported in [1 1] were obtained only in 
the case of one primary tumor. However glioma cells, being highly 
invasive, are more likely to have infiltrated into the brain tissue 
and formed multifocal or satellite lesions at time of diagnosis. In a 
recent univariate analysis of 1 5 imaging features of patients with 
GBM [61], Pope el al. found that multifocality/satellite lesions, in 
addition to the noncontrast-enhancing tumor and edema, were 
statistically significant prognostic indicators. It is therefore 
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Figure 7. Effect of Chase-ABC injections on killing infiltrating glioma cells. (A) Initial distribution of virus (v) and tumor cells (invasive cells 
outside the tumor core and proliferative cells within the ellipsoidal tumor bulk; x). (B-H) Initial profiles of Chase-ABC (upper panels) and infected cell 
density at final time t = 60h (lower panels). Viruses were injected on the periphery of the tumor bulk and Chase-ABC was injected at the different sites 
with equal strength. Number of injection sites: 1 (B), 2 (C), 3 (D), 4 (E), 5 (F), 6 (G), 7 (H). All parameters are fixed as in Table 1. 
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interesting to determine the outcome of applying viral therapy in 
the case where primary and secondary tumors are present. This is 
particularly challenging in the case where as the size of one of the 
secondary tumor is increasingly small (representing small satellite 
lesions) or when the two tumors are not too close together or not 
too far apart. In this section we consider these two situations. 

We first investigate the effect of the secondary tumor size on 
virus infection of primary and secondary tumors in the absence 
and presence of Chase-ABC. Figures 10B,C show initial distribu- 
tions of uninfected tumor cells (primary tumor with a ellipsoidal 
shape and spherical secondary tumor as well as invasive cells 
outside these tumors) when the radius of the secondary tumor (R s ) 
is relatively large (R S = 0.12; Figure 10B) and small (,R S = 0.03; 
Figure IOC). Viruses (v) were injected at five sites outside the 
primary and secondary tumors (Figure 10A). In the absence of 
Chase (Chase-ABC(— )), there are no much differences in relative 
infection rates of both primary and secondary tumors for the large 
(R S = 0A2; Figure 10E) and smaU (7^ = 0.03; Figure 10F) size of 
the secondary tumor. In fact, relative infected areas does not 
change much for various radii (R s = 0.03, 0.04, 0.05, 0.07, 0.08, 
0.09, 0.1, 0.1 1, 0.12) of the secondary tumor in both primary (left 
panel; white bar in Figure 10J) and secondary (right panel; white 



bar in Figure 10J) tumors in the Chase-ABC(— ) case. In the 
presence of Chase-ABC (Chase-ABC(+)), there are significant 
improvements relative to the Chase-ABC(— )) case in infection 
rates of both primary (left panel; black bar in Figure 10J) and 
secondary (right panel; black bar in Figure 10J) tumors for the 
fixed R s . As the secondary tumor size (R s ) is increased, the relative 
infected area of the secondary tumor is decreased but the infection 
rate of the primary tumor stays almost constant. For example, 
while there are no much differences in infected areas of the 
primary tumor at final time t = 60h for a large (_R S = 0.12; 
Figure 10H) and small (i?, = 0.03 in Figure 101) secondary tumor, 
a large portion of the larger secondary tumor was not infected 
(upper right comer in Figure 10H) and all areas of the smaller 
secondary tumor were infected (upper right comer in Figure 101). 
From the above experiments we conclude that with Chase 
treatment, the relative infected area increases as the size of the 
secondary tumor decreases. This is due to efficient spread of Chase 
throughout the relatively smaller satellite, followed up by the 
aggressive and effective degradation of CSPG, and thus resulting 
in much less physical resistance for the virus spread. 

As a second test, we investigate the effect of the distance 
between primary and secondary tumors on virus infection of these 
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Figure 8. Optimal strategies of killing infiltrating glioma cells with various Chase-ABC injections. (A) Relative infected area for various 
number of sites and injection strength of the Chase-ABC. Viruses were injected on the periphery of the tumor bulk and the Chase-ABC was injected at 
the different sites with different strength. (B) Relative infected area for invasive (white), core (gray), and total glioma cells in eight cases in Figure 7. All 
parameters are fixed as in Table 1. 
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tumors in the presence of Chase-ABC (Chase-ABC(+)). Figur- 
es 11B,C show initial distributions of uninfected tumor cells 
(primary tumor with a ellipsoidal shape and spherical secondary 
tumor as well as invasive cells outside these tumors; x) when the 
distance between the primary and the secondary tumors (d s ) is 
relatively large (d s = 0.5; Figure 1 IB) and small (<i T = 0.41 ; 
Figure 11C). Here, the distance between two tumors was defined 
as the distance between the center points of these tumors, i.e., 
d s = |x p — x s | where x p and x s are the center points of the primary 
and the secondary tumors, respectively. Viruses (v) were injected at 
five sites outside the primary and secondary tumors (Figure 1 1A). 
Initial profile of Chase- ABC (Oj is shown in Figure 1 ID. In order 
to see the effect of Chase-driven enhancement of ECM 
degradation on the efficacy of virus spread on both tumors, two 
different ECM degrading rates (X s 4l = for the control (normal) 
case; = IOA41 for the high degradation rate) for the secondary 
tumor were considered while the rate for the primary tumor was 
fixed =Xa\). Here, A41 is the control parameter in the model. 
Figure 1 1 1 shows relative infected areas of the primary and 
secondary tumors for various d s (d s = 0.5, 0.47, 0.44, 0.41). In both 
cases, relative infected areas for both primary and secondary 
tumors are increased as the distance between these tumors (d s ) is 
decreased (white bars in Figure 1 1). When the ECM-degrading 
strength of the secondary tumor is increased by 10-fold 
(a s ai = IO241), the relative infected areas of the secondary tumor 
is significantly increased regardless of distance between two tumors 
(right panel in Figure 1 II). However, the model predicts that the 
virus infection of the primary tumor in response to the larger A4] 
depends on d s in a non-linear fashion (left panel in Figure 1 II). 
While the enhanced strength of ECM-degradation of the 
secondary tumor (AJj) slightly decreases the relative infected area 
of the primary tumor for a relatively larger fixed d s 
[d s = 0.5,0.47,0.44), it promotes virus spread in the primary tumor 
for the close secondary tumor (£4 = 0.41). This is because the 
relative spread of Chase toward the primary tumor is reduced due 
to increased degradation of the secondary tumor followed by the 
increased spread speed of Chase-ABC toward the secondary 



tumor. When the secondary tumor is close to the primary tumor 
(d s = 0.41), the high viral activity in the secondary tumor may 
induce the fast viral spread to the neighboring primary tumor 
leading to increased viral infection at the primary tumor site. For 
example, one can see the broad infected areas in the secondary 
tumor and bridge-like spread of viruses between the primary and 
secondary tumors when these tumors are close each other 
(c/, = 0.41; Figure 11H) while the limited spread of viruses is 
observed in the remotely located primary tumor (d s = 0.5; 
Figure 11G). Invasive cells are more susceptible to virus infection 
due to free access of viruses and these infiltrating cells in the gap 
between the two tumors also contribute to the larger infection area 
of the primary tumor by accelerating the generation of available 
viruses in the neighboring region close to the primary tumor, 
acting as a bridge for virus transport between two tumors. Control 
cases for these remote and close satellites are shown in 
Figures 1 IE, 1 IF, respectively. From the above experiments we 
see with Chase treatment, the relative infected area increases as 
the distance from the satellite to the primary tumor decreases. This 
is due to the advantage of having larger virus population near the 
primary lesion from the active viral infection, and hence heavy 
viral duplication, on the outer rim of the primary tumor. 

Combining the above two experiments it follows that the 
relative infected area will be relatively large if the secondary tumor 
is small and close to the primary tumor; this suggests that Chase 
may perhaps be considered as anti metastatic enhancement of 
viral therapy. 

Discussion 

Genetically manipulated oncolytic viruses showed great poten- 
tial of destroying cancer cells with quick replication ability but with 
minimal damage to normal healthy cells [4]. Despite its great 
potential and popularity in China, US, and Europe [5], its 
inefficient OV spread over the entire tumor ECM has been 
recognized as a major challenge in anti-tumor efficacy [5,62]. 

In gliomas, systemic metastasis is rare and 90% of gliomas recur 
locally (less than a few cm from the resected area) [63,64], making 
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Figure 9. Effect of intravenous injection of viruses on killing infiltrating glioma cells with no Chase-ABC(— ) and fixed Chase-ABC(+) 
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gliomas good targets of virus therapy [56]. Accumulated CSPGs in 
the glioma ECM is associated with tumor growth /invasion and 
angiogenesis [29] by increasing tortuosity of the extracellular space 
and interstitial pressure within the tumor [11], resulting in 
blocking transport of large therapeutic agents [22]. Tumor 
proliferation and dispersion can be inhibited when ECM 
macromolecules are reduced [30,31]. Dmitrieva et al. [11] found 
that OV spread was significantly enhanced by OV Chase 
treatment, making ChaseABC one of important therapeutic 
molecules that can improve interstitial transport of therapeutic 
agents into the tumor. In this study, Chase expression did not 
interfere with viral cytotoxicity. Virus penetration in a local 
microenvironment may depend on many factors such as ECM 
composition and deposition of collagen and fibrinogen from highly 
permeable angiogenic vessels [23], and constant remodeling of 
ECM [65]. 

In this paper we model the phenomena of Chase-assisted virus 
spread using a system of partial differential equations. The model 
variables include tumor cells (infected or uninfected by virus), dead 
tumor cells, OV, Chondroitinase, and the ECM component. 
The model simulation are in qualitative agreement with the 



experimental results of Dmitrieva el al. [11]. It is a common 
practice to inject the virus into several different location within or 
outside the tumor mass. The mathematical model can be used to 
suggest testable hypotheses on preferred locations for such 
injections. In the absence of Chase, viral spread is very limited 
after injection of the oncolytic virus and only the cells on the 
periphery of the spheroid are infected leaving bulk of tumor cells 
uninfected at the center of tumor mass. On the other hand, viral 
penetration toward the center of tumor spheroid is improved with 
assistance of Chondroitinase Chase and much larger portion of 
tumor cells are infected near the center of the tumor spheroid. 
This improvement is due to the ability of Chase to degrade the 
thick network of ECM allowing virus to enter more easily into the 
neighboring site. This viral penetration can be improved when 
smarter virus is used, i.e., when virus is equipped with an ECM 
degrading enzyme through infected cells [51]. This effect was 
tested in our model and its improved efficacy was illustrated in our 
model simulations (Figures 2). Therefore, development of smart 
virus armed with bacterial chondroitinase ABC or other powerful 
ECM degrading enzyme may provide a pivotal way of promoting 
virus spread within the tumor by increasing enzyme activity, thus 
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killing off tumor cells. However, we note that the use of these 
degrading enzymes may cause serious complications such as 
hemorrhagic necrosis of brain, neurodegenerative diseases [26], 
and optic glioma growth [2 7] . Treatment failures in glioblastoma 
are partially attributed to cellular heterogeneity [66,67] and local 
microenvironment [36] . We also found that the microenvironment 
plays a significant role in viral spread as demonstrated in Figure 4, 
and this has implication for brain tissue with variable C-diffusion 
coefficient Dq- Our model predicts that the virus spread is quite 
selective in given microenvironment they are imposed in and 
location-specific delivery of virus is needed in order to widen 
infected area and improve the efficacy of the virus treatment. 

Infiltrating glioma cells create major challenge in glioma 
treatment due to their invisibility for standard medical devices 



such as MRI [68]. These cells may infiltrate the narrow 
intercellular space between normal cells in the brain tissue by 
deforming cell body and nucleus with effective use of myosin II 
motor [69]. Therefore, it is important to develop an efficient 
delivery method to target and kill dispersed tumor cells [58]. In 
orthotopic glioma models in rodent brain, systemic application of 
intravenous vesicular stomatitis virus (VSV) may be effective on 
selective targeting and infecting remote satellite glioma cell clusters 
as well as the tumor bulk [58,70]. We tested the efficacy of virus 
therapy on eradicating those cells in addition to killing the tumor 
cells in the main tumor mass. Under the biochemical and 
biophysical constraints, optimal virus therapy with Chase-ABC 
treatment would be also effective on killing these invasive cells as 
long as they are staying within the local areas. Our study found 
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that the injection of Chase-ABC in addition to injection of virus 
will be effective not only in infecting and killing the tumor bulk but 
also in eradicating infiltrating glioma cells even in the presence of 
the low density of blood vessels. However, this may not be effective 
on killing these migratory glioma cells in the far field even with the 
optimal strategies [46] . In addition, entry of even small biological 
agents into brain tissue is usually more complex than other organs 
because of the blood-brain barrier (BBB) [59,60,71], adding more 
obstacles in viral therapy. Some of replicating therapeutic viruses 
have not been able to disseminate to some parts of tumor [72,73]. 
Viral therapy may affect only local areas near blood vessels [73]. 
For instance, a histological analysis from transduced human 
glioma showed low transduction efficiencies with a retrovirus LacZ 
marker gene and adenoviral LacZ marker gene in some regions of 
the tumors while the neighboring areas showed much lower 



efficiency. These studies indicate that multiple injections are 
needed to infect fully the tumor bulk or the walls of the tumor 
cavity [74—76] . Some improved viruses were shown to target and 
kill cells in different regions inside the tumor via vessel leakage, i.e., 
virus seeding at multiple sites [58]. Our study predicts the 
improvement in an intravenous viral delivery when appropriate 
levels of Chase-ABC were used in an optimal way. Some viruses 
also depend on specific genetic defects such as RAS or myc 
oncogene activity, cell cycle status, p53 status, or hypoxic 
environment [77-79]. Therefore, optimized patient-specific strat- 
egies would be necessary in order to get better clinical outcomes, 
without causing catastrophic results. From our in silica experi- 
ments, we found that the size of tumor satellites and distance 
between the primary tumor and secondary ones may determine 
the efficacy of viral therapy with Chase treatment: (i) the relative 
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infected area increases as the size of the secondary tumor decreases 
because of effective degradation of CSPG and reduction in 
physical resistance of the viral spread, (ii) the relative infected area 
increases as the distance from the satellites to the primary tumor 
decreases because of the advantage of having relatively large virus 
population near the primary lesion from the active viral infection, 
and hence heavy viral duplication, on the outer rim of the primary 
tumor. 

In the current model, we did not take into account many factors 
such as cell migration [1], angiogenesis, stromal cells within the 
tumor [80], constant remodeling of ECM from blood flow, 
immune system control or growth factors [81] such as epidermal 
growth factors (EGF) [82] and fibroblast growth factors (FGF) [83] 
that may play significant roles in regulating efficacy of oncolytic 
virus therapy. In fact, those growth factors were used as a tag or 
associated molecule for targeting cancer cells by oncolytic virus 
[84]. For example, Verheije el al. [85] showed that infection of 
malignant human glioblastoma U87AEGFR cells resulted in 
release of progeny virus and efficient elimination of cancer cells in 
vitro [85] . Many research groups studied strategies for targeting of 
polymer-coated adenovirus to the EGFR [86]. How temporal- 
spatial dynamics of these growth factors, such as fluctuations, 
affect glioma cell invasion/proliferation or efficacy of viral therapy 
is largely unknown. Also, despite all advantages of oncolytic virus 
as carrier for gene therapy or themselves, the clinical application is 
limited due to practical safety issues such as immune response [87] 
and tropism of virus [88], and rather systemic administration is 
necessary [84,88]. Also, we might get better insight by considering 
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both free (interstitial) and bound virus [89] instead of one 
population of virus. However, the mathematical model developed 
herein is a first step toward incorporating these factors, which we 
hope to address in future work. We finally note that other 
approaches such as a hybrid model [55,90,91] can be useful to 
describe detailed mechanical effects such as virus penetration of a 
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from surgery due to low density, and eventually recur [1]. 
Individual cell migration might depend on many factors such as 
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